All-optical generation of states for "Encoding a qubit in an oscillator 
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Both discrete and continuous systems can be used to encode quantum information. Most quantum computation 
schemes propose encoding qubits in two-level systems, such as a two-level atom or an electron spin. Others 
exploit the use of an infinite-dimensional system, such as a harmonic oscillator. In "Encoding a qubit in 
an oscillator" [Phys. Rev. A 64 012310 (2001)], Gottesman, Kitaev, and Preskill (GKP) combined these 
approaches when they proposed a fault-tolerant quantum computation scheme in which a qubit is encoded in 
the continuous position and momentum degrees of freedom of an oscillator. One advantage of this scheme is 
that it can be performed by use of relatively simple linear optical devices, squeezing, and homodyne detection. 
However, we lack a practical method to prepare the initial GKP states. Here we propose the generation of an 
approximate GKP state by using superpositions of optical coherent states (sometimes called "Schrodinger cat 
states"), squeezing, linear optical devices, and homodyne detection. 



The Gottesman, Kitaev, and Preskill (GKP) scheme 
[1], constitutes a type of linear optical quantum com- 
puter as are other schemes based on the proposals of 
Knill, Laflamme, and Milburn [2] and schemes based on 
the proposal of Ralph et al. [3] . In the GKP scheme, the 
qubit is encoded in the continuous Hilbert space of an 
oscillator's position and momentum variables. The GKP 
scheme is applicable to any type quantum harmonic os- 
cillator, but we focus on an optical implementation in 
traveling modes. The GKP encoding provides a natural 
error-correction scheme to correct errors due to small 
shifts (applications of the displacement operator) on the 
conjugate quadrature variables x and p [1]. 

The ideal GKP logical qubit state, |0), is defined 
as a state whose cc-quadrature wave function is an infi- 
nite series of delta-function peaks occurring whenever 
x = 2y / 7ts for all integers s, while the ideal |1) x- 
quadraturc wave function is displaced a distance ^/ir 
from the |0) state. Since these states arc unphysical, 
GKP described approximate states whose ^-quadrature 
wave function is a series of Gaussian peaks with width 
A, contained in a larger Gaussian envelope of width 1/fc. 
The approximation of |0) has the wave function 
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where N is a normalization factor. If A and k arc small, 
then ipGKP will better approximate an ideal GKP state, 
and the wave function will have many sharp Gaussian 
peaks contained in a wide envelope. We can think of 
the deviation from an ideal GKP state as corresponding 
to nonzero probability that the state has suffered from 
errors causing displacement in the x or p variables. If all 
displacements are smaller than y/ir/6, then the errors will 
not increase during the error correction protocol. GKP 
states with A < 0.15 and k < 0.15 will have a probability 
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Fig. 1. Approximate GKP state's x-quadrature wave 
function. This shows the logical state ipGKp{x) with 
A = k = 0.15. 



greater than 0.99 of being free of shift errors larger than 
\/tt/6 [4]. Fig. 1 shows an example of an approximate 
GKP state's x-quadrature wave function. 

Preparing GKP states is a difficult task, and no exper- 
iment has yet demonstrated preparation of such states. 
GKP proposed preparing these states by coupling an op- 
tical mode to an oscillating mirror in [1]. Another pro- 
posal was made by Travaglione and Milburn in [5] , where 
the qubit states are prepared in the oscillatory motion 
of a trapped ion rather than the photons in an optical 
mode. Pirandola et al. [6] discusses the preparation of 
optical GKP states by use of a two mode Kerr inter- 
action followed by a homodyne measurement of one of 
the modes. The same authors describe two proposals for 
generating GKP states in the position and momentum 
of an atom by using a cavity mediated interaction with 
light in [7,8]. 

Here we propose the generation of an approximate 
GKP state by using superpositions of optical coherent 
states ("cat states"), linear optical devices, squeezing, 
and homodyne detection. The basic idea is: first, pre- 
pare two cat states (each of which contains two Gaussian 
peaks in its x-quadrature wave function), squeeze both 
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cats (to reduce the width of the Gaussian peaks), inter- 
fere them at a beam splitter, then perform homodyne 
detection on one of the beam splitter's output ports. 
Depending on the measurement result, we will find an 
approximate GKP state (with three Gaussian peaks) in 
the beam splitter's other output port. This procedure 
can be repeated to produce states with larger numbers 
of Gaussian peaks. 

A cat state is a superposition of coherent states such 
as: 

I / i w I ~ g) + I") to\ 
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where a is the amplitude of the coherent state, which 
may be complex, but we assume it is real below. Several 
experimental proposals to create cat states are reviewed 
in [9] . Cat states of this form have been created in several 
experiments with \a\ up to 1.75 and fidelities of 0.6 to 
0.7 [10-15]. 

Cat states' ^-quadrature wave functions are superposi- 
tions of Gaussian peaks. To simplify notation, we denote 
a Gaussian with 
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These Gaussians represent wave functions, not probabil- 
ity distributions, so the unnormalized vacuum state is 
G(x, 1,0). Suppose two modes, labeled 1 and 2, contain 
states with unnormalized wave functions G{x\ 1 V, n\) 
and G(x2, V,fj, 2 ). Modes 1 and 2 meet at a beam splitter 
with transmissivity 1/2, which performs the transforma- 
tion: 
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After the beam splitter, we use a homodyne detector 
to measure mode 2's p-quadrature. In the case that the 
measurement result is p 2 = 0, this entire procedure pro- 
duces the transformation 

G(x 1 ,V,m)G(x 2 ,V,fi 2 ) WG(x u V 7 ^J^). (5) 

We can write Eq. (2) in the x-quadrature basis as: 

=G(x,l,-V2a) + G{x,l,V2a), (6) 

where we use the tilde to signal that the state is not 
normalized. We now squeeze this state by an amount £, 
obtaining 

■4> sq cat(x,a, C) = G(x, e~ 2C , -V^ae^) 

+G(x,e- 2C ,V2ae- c ). (7) 

We will choose the cat state amplitude to be a = 
y/ne^, where m is an integer greater or equal to 1, 



which we will later use to count iterations of our scheme. 
With these choices, 
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Suppose we have two copies of this squeezed cat in modes 
1 and 2. They meet at a beam splitter with transmissivity 
1/2, and the p-quadrature of mode 2 is measured to be 
p 2 = 0. If we choose m = 1, the resulting unnormalized 
state of mode 1 is 
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which we call the first binomial state and it is similar 
to an approximate GKP logical qubit 0, except only the 
central three peaks are present. The width of those peaks 
is controlled by the amount of squeezing £ applied to the 
original cat states. 

Consider the order to binomial state given by 
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This state is a series of Gaussian peaks separated by 2y / 7r 
along the x-quadrature axis, and the amplitudes of the 
peaks are given by the (2 m ) th row of Pascal's Triangle 
(where row contains only 1). We will show that given 
two copies of the order to binomial state, one can make 
the m + 1 order binomial state. We begin with the state 
of modes 1 and 2: 
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These two modes meet in a beam splitter of transmissiv- 
ity 1/2, and we measure the p-quadrature, obtaining the 
result p 2 = 0. The new state is given by applying Eq. (5) 
to $(xi, C, m) (3{x2,C, 1 to). The result is 
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After a little algebra and application of Vandcrmont's 
identity, we obtain 

E ( q )G(x l7 e-^,V2^(q-2 m ))- (13) 

This state is equivalent to /3(xi,£, m + l)i except that 
the Gaussian peaks are separated by only V^^/n rather 
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Fig. 2. Left: Probability for measuring p 2 ~ R a.s func- 
tion of R. Right: ^-quadrature wave function for the state 
(3{x, £, 1). In both cases a = y^e^, and Q = 1.9 is chosen 
to produce Gaussian peaks with the same width as shown 
in Fig. 1. This ( is equivalent to -16 dB of quadrature 
noise power reduction in conventional squeezing experi- 
ments. 




Fig. 3. Wave function of binomial state j3(x, 3), a closer 
approximation of the logical GKP state. Again here 
C = 1.9. Creating this state would require at least eight 
cat states, each with a = 2 y /ire'>. The /3(x, £, 3) state also 
has small peaks at ±8y/n, which are not visible here. 

than 2 v / 7r. We can compensate for this shrinking of the 
spacing between the Gaussian peaks, if we begin the pro- 
cedure with two binomial states with spacing of 2-s/Vtt. 
As m increases, these binomial states will approach the 
shape of a series of Gaussian peaks in a Gaussian enve- 
lope, like ■4> GKP {x). 

In Fig. 2, we plot the probability for measuring a cer- 
tain value P2 = R as a function of R when making the 
m = 1 binomial state from two cat states. If we mea- 
sure R — 0, we obtain the wave function as shown in 
Fig. 2, whose Gaussian peaks' widths are determined by 
the degree of squeezing applied to the initial cat states, 
and whose heights arc proportional to the second row of 
Pascal's triangle (1, 2, 1), as given by Eq. (10). In Fig. 3 
we show the m = 3 binomial state. 

Creating the order m binomial state requires a min- 
imum of 2 m cats, but the true number may be much 
larger, because we require that p 2 = at each measure- 
ment event. However, initial investigations indicate that 
some cases in which p 2 can be recovered by applying 
a p-quadrature displacement whose size depends on the 
measurement result. We plan to explore this further in 
a future work. Creating high quality states will require 
larger squeezing £ and higher order m binomial states. 
We also plan to investigate the performance of these bi- 
nomial states in the GKP encoding scheme as a function 



of £ and m. 

In this paper we began construction of GKP with a 
source of cat states. The most popular method to make 
cat states is by subtracting photon(s) from a squeezed 
vacuum state. It may be possible to alter the photon 
subtraction scheme to benefit our method to make GKP 
states. 

Although our scheme is built of apparently simple, 
well understood optical operations, it will be difficult to 
achieve in an experiment. We expect that matching the 
transverse and longitudinal shapes of all of the optical 
modes, especially during the squeezing stage [16,17], will 
be very difficult. 

We thank Adam Meier, Yanbao Zhang, and Emanuel 
Knill for helpful discussion and comments. H. M. Vas- 
concelos thanks the PNPD program (CAPES) and FUN- 
CAP for financial support. L. Sanz thanks INCT-IQ for 
financial support. This paper is a contribution by the 
National Institute of Standards and Technology of the 
United States of America and not subject to USA copy- 
right. 

References 

1. D. Gottesman, A. Kitaev, J. Preskill, "Encoding a qubit 
in an oscillator", Phys. Rev. A 64, 012310 (2001). 
http://arxiv.org/abs/quant-ph/0008040. 

2. E. Knill, R. Laflamme, G. Milburn, "A scheme for effi- 
cient quantum computation with linear optics," Nature 
409, 46-52 (2001). 

3. T. Ralph, A. Gilchrist, G. Milburn, W. Munro, 
S. Glancy, "Quantum computation with optical co- 
herent states," Phys. Rev. A 68, 042319 (2003). 
http://arxiv.org/abs/quant-ph/0306004. 

4. S. Glancy, E. Knill, "Error analysis for encoding a 
qubit in an oscillator," Phys. Rev. A 73, 012325 (2006). 
http://arxiv.org/abs/quant-ph/0510107. 

5. B. C. Travaglione, G. J. Milburn, "Preparing encoded 
states in an oscillator," Phys. Rev. A 66, 052322 (2002). 
http://arxiv.org/abs/quant-ph/0205114. 

6. S. Pirandola, S. Mancini, D. Vitali, P. Tombesi, "Con- 
structing finite-dimensional codes with optical contin- 
uous variables," Europhys. Lett. 68, 323-329 (2004). 
http://arxiv.org/abs/quant-ph/0402202. 

7. S. Pirandola, S. Mancini, D. Vitali, P. Tombesi, "Gen- 
erating continuous variable quantum codewords in the 
near-field atomic lithography," J. Phys. B.39, 997-1009 
(2006). http://arxiv.org/abs/quant-ph/0510053. 

8. S. Pirandola, S. Mancini, D. Vitali, P. Tombesi, 
"Continuous variable encoding by ponderomotive in- 
teraction," Eur. Phys. J. D 37, 283-290 (2006). 
http: / /arxiv.org/abs / quant-ph/0503003 

9. S. Glancy, H. M. Vasconcelos, "Methods for 
producing optical coherent state superposi- 
tions," J. Opt. Soc. Am. B 25 712-733 (2008). 
http://arxiv.org/abs/0705.2045 

10. A. Ourjoumtsev, R. Tualle-Brouri, J. Laurat,P. Grang- 
ier, "Generating Optical Schrodinger Kittens for Quan- 
tum Information Processing," Science 312, 83-86 (2006). 

11. J. S. Neergaard-Nielsen, B. Melholt Nielsen, C. Hettich, 
K. M0lmer, E. S. Polzik, "Generation of a Superposi- 



3 



tion of Odd Photon Number States for Quantum Infor- 
mation Networks," Phys. Rev. Lett. 97, 083604 (2006). 
http://arxiv.org/abs/quant-ph/0602198. 

12. K. Wakui, H. Takahashi, A. Furusawa, M. Sasaki, "Pho- 
ton subtracted squeezed states generated with period- 
ically poled KTiOP0 4 ," Opt. Express 15, 3568-3571 
(2007). http://arxiv.org/abs/quant-ph/0609153. 

13. A. Ourjoumtsev, H. Jeong, R. Tualle-Brouri, P. Grang- 
ier, "Generation of optical 'Schrodinger cats' from pho- 
ton number states," Nature 448, 784-786 (2007). 

14. H. Takahashi, K. Wakui, S. Suzuki, M. Takeoka, K. 
Hayasaka, A. Furusawa, and M. Sasaki, "Generation 
of Large-Amplitude Coherent-State Superposition via 
Ancilla-Assisted Photon Subtraction," Phys. Rev. Lett. 
101, 233605 (2008). http://arxiv.org/abs/0806.2965. 

15. T. Gerrits, S. Glancy, T. S. Clement, B. Calkins, A. E. 
Lita, A. J. Miller, A. L. Migdall, S. W. Nam, R. P. Mirin, 
E. Knill, "Generation of optical Schrodinger cat states 
by number-resolved photon subtraction from squeezed 
vacuum," http://arxiv.org/abs/1004.2727. 

16. W. Wasilewski, A. Lvovsky, K. Banaszek, C. Radzewicz, 
"Pulsed squeezed light: Simultaneous squeezing of 
multiple modes," Phys. Rev. A 73, 063819 (2006). 
http://arxiv.org/abs/quant-ph/0512215. 

17. O. Cohen, J. S. Lundeen, B. J. Smith, G. Puentes, P. J. 
Mosley, I. A. Walmsley, "Tailored Photon-Pair Gener- 
ation in Optical Fibers," Phys. Rev. Lett. 102, 123603 
(2009) . http: / / arxiv.org/abs/0809.0071v2. 



4 



